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ABSTRACT 

We describe a search of archival data from the Burst and Transient Source Experi- 
ment (BATSE). The purpose of the search is to find astronomically interesting transients 
that did not activate the burst detection (or "trigger") system onboard the spacecraft. 
Our search is sensitive to events with peak fluxes (on the 1.024 s time scale) that are 
lower by a factor of ~ 2 than can be detected with the onboard burst trigger. In a 
search of 345 days of archival data, we detected 91 events in the 50-300 keV range that 
resemble classical gamma ray bursts but that did not activate the onboard burst trigger. 
We also detected 110 low-energy (25-50 keV) events of unknown origin which may in- 
clude activity from SGR 1806—20 and bursts and flares from X-ray binaries. This paper 
gives the occurrence times, estimated source directions, durations, peak fluxes, and flu- 
ences for the 91 gamma ray burst candidates. The direction and intensity distributions 
of these bursts imply that the biases inherent in the onboard trigger mechanism have 
not significantly affected the completeness of the published BATSE gamma ray burst 
catalogs. 

Subject headings: gamma rays: bursts 
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1. Introduction 

Since 1991 April 19 the Burst and Transient Source 
Experiment (BATSE) on the Compton Gamma Ray 
Observatory (CGRO) has been detecting gamma ray 



unprecedented sensitivity ( 


Fishman et al. 198E ; 


Fish- 


man ct al. 1994b 




Meegan et al. 1996 


). The 1122 



GRBs in the 3B catalog show an isotropic angu- 
lar distribution and a spatially inhomogeneous inten- 



sity distribution (Meegan ct al. 1992; Meegan et al 



1996 ). Despite extensive analysis, however, the ori- 
gin of GRBs remains unknown; see recent reviews by 
Fishman & Meegan (1995), Briggs (1995), and Hart- 
mann (1995). 



L he detection ot Gttbs and other high-energy tran- 



sients with tiAl "SE is controlled by a real-time burst 
detection algorithm running onboard the spacecraft 
( Fishman ct al. 1989| ). The onboard computer con- 
tinuously monitors the count rates in each of the eight 
Large Area Detectors (LADs) . When the count rates 
exceed a certain threshold, the computer signals a 
"burst trigger" and data are collected at high tem- 
poral and spectral resolution for a limited time in- 
terval. Even in the absence of a burst trigger, how- 
ever, data are recorded at lower resolution in the con- 
tinuous data types. For most of the mission, the 
criteria for a burst trigger have been that the 50- 
300 keV count rates in two detectors simultaneously 
increase by more than 5.5 times the expected root- 
mean-square background fluctuations on any of three 
time scales: 64 ms, 256 ms, or 1024 ms. The aver- 
age backgroun d rate for each detect or is recomputed 
every 17.408 s ( [Fishman et al. 1989[ ). 

By definition all of the GRBs listed in the IB, 
2B, and 3B c atalogs satisfy the requirements for a 
burst trigger (Fishman et al. 1994b; Meegan et al 



1996). Other transient phenomena that are unrelated 
to GRBs can also lead to a burst trigger. Examples 
include solar flares, terrestrial magnetospheric distur- 
bances, bursts and flares from X-ray binaries, and ac- 
tivity from soft gamma ray repeaters (SGRs). Such 
events are classified appropriately by the BATSE 
team. 

A GRB or other transient phenomenon may have 
characteristics such that it does not lead to a burst 
trigger onboard the spacecraft but it nevertheless 
leaves a statistically significant signal in the continu- 
ous data. For example, a GRB or other transient may 
be too faint to achieve the necessary statistical signif- 



icance for a trigger; it may have a time profile that bi- 
ases the onboard background average; or it may have 
too few counts in the 50-300 keV range. 

A burst can also occur while the onboard trigger 
is disabled for technical reasons. Following a burst 
trigger, the high resolution data collected during the 
burst accumulation interval are gradually telemetered 
to the ground during the following 90 minutes. Dur- 
ing this read out period the onboard burst trigger is 
disabled on the 256 ms and 1024 ms time scales, and 
the 64 ms threshold is set to the maximum rate of 
the burst being read out. The onboard burst trigger 
is also disabled when the spacecraft passes through 
regions with a high probability of triggering on atmo- 



spheri c particle precipitation events (Fishman et al 



1994b[) 



In this paper, we describe a retrospective search 
of the archival continuous data from BATSE for sta- 
tistically significant GRBs and other transients that 
did not cause a burst trigger onboard the space- 
craft. A search for these "non-triggered" (or "untrig- 
gered") events in the 50-300 keV range is expected 
to find GRBs that are generally fainter than those 
cataloged previously. A concurrent search for non- 
triggered events in the lowest discriminator channel 
(25-50 keV) is expected to find activity from other 
astronomical sources, including bursts and flares from 
X-ray binaries and activity from SGRs. This ongoing 
project is an extension of previous work by Rubin et 
al. (1993), Van Paradijs et al. (1993), and Kommers 
et al. (1996). Other retrospective searches for GRBs 
in the BATSE data (using techniques different from 
those described here) have been discussed by Skelton 
& Mahoney (1994) and by Young et al. (1996). 

2. Search Algorithm 

A retrospective search of archival data can take ad- 
vantage of burst detection algorithms that would have 
been impractical to implement onboard the space- 
craft. The choice of a detection scheme to look for 
non-triggered events therefore involves trade-offs be- 
tween the detection efficiency of the method for a 
given class of events and the resources (both com- 
putational and human) needed to implement it. In 
this section we describe a search algorithm which is 
loosely based on the one used onboard the spacecraft 
but which has proved more sensitive. The next sec- 
tion will discuss its efficiency for detecting transients 
with certain characteristics. 
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We will refer to all events detected by our off- 
line search of archival data as "laboratory triggers." 
The events previously detected by the onboard burst 
trigger mechanism will be called "onboard triggers." 
Some onboard triggers will be flagged by our off-line 
search and so they will also be laboratory triggers. 
Events that were detected only by our off-line search 
will be called "non-triggered events." 

For most of the mission the onboard trigger crite- 
rion has required that the count rate in two detectors 
simultaneously increase by at least 5.5<tb above the 
nominal background level, where o~b is the standard 
deviation of the expected background counts due to 
counting statistics. As a result, the BATSE detectors 
provid e anisotropic sky exposure over short time pe - 
riods (|Fishman et al. 1989| ; |Fishman ct al. 1994b| ). 



The cosine-like change in the detectors' effective area 
with source viewing direction causes the onboard trig- 
ger to be less sensitive to faint bursts with directions 
directly in front of one of the detectors than to ones 
with directions mid-way between two detector nor- 
mals ( Brock et al. 1991 ) . For example, a IOctb event 
occurring directly in front of a detector may produce 
only a 3.5(7b signal in the second most brightly il- 
luminated detector and it would thus fail to trigger 
onboard. On the other hand, the same event inci- 
dent along a direction mid-way between two detector 
normals would register approximately 7. las in both 
detectors and would comfortably cause an onboard 
trigger. 

The onboard trigger mechanism relies on a back- 
ground rate that is computed during a 17.408 s time 
interval occurring before the time bin being tested. 
A rising or falling background can therefore bias the 
background estimate to be too low or too high, re- 
spectively. A slowly rising transient may itself bias 
the background estimate upwards to such an extent 
that it fails to cause an onboard burst trigger, even 
though it is otherwise intense enough to be above the 
minimum detection threshold. 

Our retrospective search procedure partially com- 
bats the directional detection anisotropy and the ris- 
ing/falling background bias. We form a time series to 
be searched by combining the relevant energy chan- 
nels from the DISCLA data. This data type provides 
the count rates in each detector integrated over 1024 
ms time bins. Four discriminator energy channels 
(numbered 1 through 4) are available: 25-50 keV, 
50-100 keV, 100-300 keV, and > 300 keV. When 
searching for GRBs, the sum of channels 2 and 3 



(50-300 keV) gives optimal sensitivity. When search- 
ing for low-energy transients such as bursts and flares 
from X-ray binaries or SGRs the lowest energy chan- 
nel (channel 1, 25-50 keV) is most sensitive. A sum 
of channels 1,2, and 3 provides a further "catch-all" 
search. After summing the appropriate energy chan- 
nels, the time series are rebinned in time (if necessary) 
to search on time scales longer than the 1024 ms DIS- 
CLA sampling period. The resulting time series are 
then searched sequentially to see if any data meet our 
laboratory trigger criteria. 

To signal a laboratory trigger, we first determine 
from the time series being searched a nominal back- 
ground level for each detector. To estimate B^ik), the 
number of background counts expected in detector d 
for the time bin k, we use a linear fit to the data in 
time bins k- N b , . . . , fc- 1 and k + l + N b , . . .,k + 2N b . 
Here the number of background bins, N b , is specified 
for each time scale. Unlike the onboard trigger, the 
off-line search uses data before and after the time bin 
being testedjj] This method reduces the bias (dis- 
cussed above) caused by slowly rising or falling back- 
ground levels. 

Temporally contiguous data are not always avail- 
able because of telemetry gaps and spacecraft pas- 
sages through regions of high particle flux, such as 
the South Atlantic Anomaly (SAA). In such cases 
the background estimate discussed here cannot be 
formed. We only test bins for which we can esti- 
mate the background with the above procedure, so 
our search is not sensitive to bursts that occur during 
the iVf, bins after, or the 2N b bins before, a data gap. 

Let Cd{k) be the measured number of counts in 
time bin k for detector d. We define the "signif- 
icance" of that detector to be Sd(k) — [Cd(k) — 
Bd{k)]/ y/ Bd{k). Our laboratory trigger criteria are 
that the two greatest values of the 8 significances 
Sd{k), call them s\ and S2, must be such that s\ > 
S2 > M and si + S2 > S. These criteria ensure that at 
least two detectors simultaneously experience a sta- 
tistically significant upward fluctuation, but they are 
also more sensitive to events incident along detector 
normals than the onboard criteria. 



lr The bin, fc, being tested is not centered in the gap of Ni, bins 
which separates the fitted intervals. This somewhat arbitrary 
choice evolved out of various triggering schemes which were 
tried, including one where the background was always esti- 
mated based on the Nf, bins immediately prior to the one being 
tested. The scheme chosen performed well on 14 days of data 
that were used to test various laboratory trigger criteria. 



3 



The values chosen for M, X, and N(, for the 
searches reported here are shown in Table |[ They 
were chosen to keep the actual number of detections 
per day of data searched (due to the activity of real 
sources) at a manageable level of about 20 per day. 

A more sensitive search could be conducted using 
other laboratory trigger criteria. For example, the 
effective detector area achieved by adding the rates 
in each set of 3 contiguous detectors would produce 
time series with higher signal-to-noise ratios for ac- 
tivity from real sources. The caveat to a more sensi- 
tive search is the corresponding increase in the rate of 
"false triggers" due to solar activity, variability from 
X-ray binaries, particle precipitation events, Earth 
occultation steps, and phosphorescence spikes (which 
occur when high-energy particles interact in the de- 
tectors). During outbursts of bright X-ray binaries, 
our laboratory trigger has detected hundreds of events 
per day due to variability from Vela X-l, A0535+26, 
and GRO J0422+32. 

For our laboratory trigger criteria, the number of 
detections expected from statistical fluctuations alone 
can be estimated by considering the "phase space" de- 
fined by S\ and S2- The number of counts in each time 
bin of the DISCLA data type is large enough that the 
Poisson statistics can be treated in the Gaussian ap- 
proximation, so that the S ( i(k) are independent and 
normally distributed with zero mean and unit vari- 
ance. The laboratory trigger criteria define a region 
A in the (si,sa) plane in which measured values s| 
and S2 will be flagged by our search. This area is 
sketched in Figure H. 

The probability P s t a t that a randomly selected set 
of 8 significances Sd will meet or exceed the laboratory 
trigger criteria is estimated by first integrating the 
bivariate normal distribution (centered on the origin) 
over the allowed area A in the (si,S2) plane. Next, 
we multiply by 8 ways to select the most significant 
detector (s%) and the remaining 7 ways to select the 
second most significant (S2). The final expression is 



2tt 



exp 



dsi ds2 



x 7 x 



(1) 

If there are N time bins searched per day, the number 
of detections expected from purely statistical fluctu- 
ations is Pstat x N per day. This quantity is shown in 
the last column of Table [l] for a representative value 
of N = 5.0 x 10 4 . Although the searches are not 
statistically independent, we would expect no more 



than 8 detections due solely to statistical fluctuations 
in a search of 100 days of data. Furthermore only a 
fraction of these will have properties consistent with 
astronomical source activity. 

In practice we find many more laboratory triggers 
than can be expected from purely statistical fluctua- 
tions. This is because the time series we are searching 
are not Poissonian. They are dominated by the ac- 
tivity of real sources. Astronomical objects, the sun, 
terrestrial photon sources, and the interaction of par- 
ticles in the detectors contribute to the count rates. 
Figure^ shows an integral distribution of Sd{k) for a 
single detector (d = 0) taken from a search of one day 
of data on the 1.024 s time scale. For comparison, 
the expectation from a normal distribution with zero 
mean and unit variance is also shown. The excess 
over what is expected from a normal distribution re- 
flects both deficiencies in the background estimation 
and the activity of real sources, although no labora- 
tory triggers were detected by our off-line search in 
this detector on the day shown. 

3. Sensitivity 

The search strategy described above is expected to 
detect bursts that were fainter than those detected 
by the onboard trigger mechanism. In this section, 
we estimate the efficiency of the search algorithm for 
detecting events with certain physical characteristics. 
In general the ability of our search strategy to de- 
tect an event depends on its peak flux (as measured 
on the search time scale), its time profile, the back- 
ground levels in the detectors, and the spacecraft ori- 
entation with respect to the source direction. Here 
we will estimate the trigger efficiency and sky expo- 
sure of our off-line search strategy assuming an event 
profile for which a single time bin completely deter- 
mines our ability to detect the event. The effects of 
a more complicated time profile are then considered 
separately for a simplified case of "slow-rising" events 
that would have biased the onboard background esti- 
mate. 

3.1. Trigger Efficiency 

We define the trigger efficiency E(P,i/,a,S,t) to 
be the probability that an event with given physical 
characteristics will satisfy the laboratory trigger cri- 
teria. The event is modeled in terms of its peak flux 
on the time scale of the search (P), power-law pho- 
ton spectral index (y), source direction in equatorial 
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coordinates (a, 6), and time of occurrence during the 
mission (t). The background rates in the detectors, 
the spacecraft orientation, and the geographic posi- 
tion of the spacecraft must also be known to estimate 
E, but these quantities are known from the data once 
t is specified. The time profile of the event is consid- 
ered to be a square pulse that occupies a single time 
bin. 

For a particular spacecraft orientation and posi- 
tion, the mean count rates expected from an event 
with given values of P, i>, a, and <5 can be com- 
puted by adding the "direct" count rates found from 
the BATSE instrument response matrix to the "scat- 
tered" count rates expected from the scattering of in- 
cident photons by the Earth's atmosphere. The in- 
strument response matrices and atmospheric scatter- 
ing rnod^laxedescribcdfurthcri^ 



(19 p5) and Meegan ct al. (1996).p| The resulting total 



count rates are multiplied by the time binning inter- 
val to obtain the expected counts (above background) 
in a single bin of the time series for each detector; we 
denote these quantities by C|. The expected back- 
ground counts, Bd, are estimated from the measured 
rates at time t and the expected significances in the 
detectors are calculated as S% = C* l /\/B~d. Let s\ and 
S2 be the the greatest and second greatest of the S% , 
respectively. 

The trigger efficiency, E, is the probability that 
a measurement of a pair (s±,S2) will meet the lab- 
oratory trigger criteria. The counting statistics im- 
ply that our measurement, (s±,S2), is drawn from 
a bivariate normal distribution with unit variances 
centered on the expected significances, (s^s^)- The 
probability that we detect the given event is there- 
fore estimated by integrating this distribution over 
the area, A (shown in Figure |l|) , in which the trigger 
criteria are satisfied: 



v = 2.0. The 9 peak fluxes were chosen to span an in- 
tensity range where the efficiency varies significantly. 
The 252 source directions are nearly isotropically dis- 
tributed on the unit sphere (Tegmark 199£). The 



4992 times per orbital precession period were chosen 
to thoroughly sample the range of background varia- 
tions. For points where the source direction is behind 
the Earth or no searchable data are available, E is 
set to zero. This calculation must be repeated for 
each time scale and energy channel combination in 
our search. 

Figure || shows £ as a function of P for events 
searched on the 1024 ms and 4096 ms time scales 
in the 50-300 keV range. The source directions and 
times of occurrence have been averaged. For compar- 
ison, the trigger efficiency on the 1024 ms t ime scale 
from the BA TSE IB catalog is also shown ( Fishman 



ct al. 1994b|) 



E 



2^ 6XP 



(si-*i) 2 



exp 



(S2 



dsi d,S2- 



(2) 

We evaluate E on & three-dimensional grid composed 
of 9 peak fluxes, 252 source directions, and 4992 times 
per orbital precession period of CGRO (i.e., every 
15 minutes). The photon spectral index is fixed at 



2 The detector response and atmospheric scattering matrices 
used in this work were provided by the BATSE team. We also 
made use of some elements of the BACODINE burst location 
code (Scott Barthelmy, private communication) which is based 
on an early version of the BATSE LOCBURST code. 



For the 1024 ms time scale, Figure y shows that our 
search is nearly complete near the BATSE threshold 
(~ 0.2 ph cm~ 2 s _1 ). Our off-line search should de- 
tect about 50% of events with peak fluxes lower by 
a factor of ~ 2 than the onboard 50% completeness 
limit. For events that maintain their peak flux for at 
least 4 s or 8 s, the searches on the 4096 ms and 8192 
ms time scales can reach even lower peak fluxes. The 
values of E(P) shown in Figure || vary by only a few 
percent between spacecraft orbital precession periods 
(52 days for CGRO). 

The sky coverage of our search is determined by 
the angular distribution of E{P, a, 5), where P and 
v are fixed and t has been averaged. If T is the to- 
tal time period covered by the data searched, then 
T x E(P, v, a, 6) gives the total amount of time that 
our search was sensitive to an event with the given in- 
tensity, spectral index, and source direction. In prac- 
tice our search covers many orbits, so the dependence 
on the equatorial right ascension (a) averages out. 
Figure ^ shows the sky exposure as a function of dec- 
lination (S) for events with P > 0.5 ph s^ 1 cm -2 on 
the 1024 ms time scale and v = 2.0. The decreased 
exposure near the celestial equator is due to Earth 
blockage. The Southern Hemisphere gets less expo- 
sure than the Northern Hemisphere due to passages 
through the SAA. 

3.2. Sensitivity to "slow-risers" 

The chances that our search detects a given event 
also depend on its time profile in our time series. For 
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example, an event with a peak flux near the detection 
threshold that is shaped like a square pulse occupying 
N time bins will have N statistical chances to meet or 
exceed the trigger criteria. In such a case the trigger 
efficiency calculated above will be an underestimate. 

On the other hand, the event profile may be such 
that it biases the nominal background rate used by 
the search algorithm and artificially raises the detec- 
tion threshold. Then the trigger efficiency calculated 
above is an overestimate. This case is more serious as 
it implies that the search could miss a population of 
such events even though their peak flux is well above 
the nominal minimum detection threshold. Both the 
onboard trigger mechanism and our off-line search are 
subject to this limitation. 

The most problematic time profile is one which 
both rises and decays slowly on a time scale long com- 
pared to the background averaging with no significant 
rapid variability. This case is difficult to distinguish 
from background variations arising from the space- 
craft environment. Neither our search nor the on- 
board trigger has appreciable sensitivity to events of 
this type. They are even unlikely to be evident in a 
close visual inspection of the count rates. 

An event that rises slowly and then either falls 
off quickly or subsequently goes into a more compli- 
cated profile can usually be identified as a transient, 
however. We will call such an event a "slow-riser." 
The onboard trigger mechanism has particular trou- 
ble with slow-risers because it can only base its back- 
ground estimate on count rates measured over some 
17.408 s interval that occurred some time during the 
34.8 s before the time bin being tested. Lingenfel- 
ter & Higdon (1996) have discussed this effect for the 
onboard burst trigger. Here we estimate the extent 
to which our off-line search algorithm is sensitive to 
slow-risers. 

The profile of an idealized slow-riser is shown in 
Figure ^[ The event is characterized by a peak ampli- 
tude, s, measured in sigmas above background, and 
a slope, r, measured in sigmas per second, where one 
sigma is one standard deviation of the expected back- 
ground fluctuations. The total duration of the event 
is s/r. To estimate our sensitivity to such events, 
we generated a grid of peak amplitudes and a grid of 
slopes. For each peak amplitude and slope, we gen- 
erated 5000 events with a background level of 1000 
counts s _1 and Poisson noise to mimic the counting 
statistics. We then used our detection algorithm to 
find the fraction of events that met our off-line trig- 



ger threshold and the fraction of events that met the 
onboard trigger threshold. 

The results are shown in Figure |^, which shows 
contours of detection probability for (idealized) events 
with maximum significance, s, and slope, r, when 
searching on the 1024 ms time scale. Evidently our 
laboratory detection algorithm is more sensitive to 
slow-risers than the onboard trigger mechanism. This 
is due both to the lower detection threshold and to the 
use of data before and after the time bin being tested 
when forming the background estimate. Slow-risers 
with no subsequent variability that are both longer 
than about 30 s and fainter than 5a (about 0.3-0.5 
photons cm -2 s _1 in the 50-300 keV range) are un- 
likely to be found by our 1024 ms off-line search. The 
4096 ms and 8192 ms searches are sensitive to longer 
events, however. 

The idealized profile used for these estimates may 
not be representative of the faintest or longest tran- 
sients since the typical profiles of those events are not 
known a priori. The simple case presented here shows 
that our search algorithm is more sensitive to some 
slow- rising bursts. 

4. Classification and Analysis of Events 

We visually inspect each laboratory trigger flagged 
by our off-line search to separate the events into useful 
categories. 

We first examine plots of the count rates in each of 
the eight detectors at the time of the laboratory trig- 
ger. This information is adequate to exclude from fur- 
ther analysis the majority of features in the data that 
are not interesting in the context of our search. Com- 
mon examples of such features are occultation steps, 
phosphorescence spikes, and magnetosphcric particle 
precipitation events. 

Occultation steps due to bright sources (such as 
Cyg X-l) rising above (or setting below) the Earth's 
limb typically appear as sustained increases (or de- 
creases) in the count rates in two or more detectors. 
This signature is not generally consistent with the 
burst-like transients we seek. 

Phosphorescence spikes due to the interaction of 
high-energy particles in the detectors typically appear 
in the DISCLA data as short (< 1.024 s), intense (> 
1000 counts s _1 ) spikes in the lowest energy channel 
(25-50 keV) of a single detector. The intense signal 
in one detector with no coincident signal in a second 
detector is inconsistent with a point source of pho- 
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tons. 

Magnetospheric particle precipitation events occur 
when particles (mostly electrons with energies of tens 
of eV to tens of MeV) that are usually trapped in 
the radiation belts by the Earth's magnetosphere are 



released into the upper atmosphere (Burgess & Inan 
19£[3j) . The LADs detect the bremsstrahlung gener- 
ated as precipitating particles interact in the atmo- 
sphere or in the spacecraft (Horack et al. 1992 and 
references therein). These events can appear in the 
data in three different ways as discussed by Horack 
et al. (1992). Events of the first kind show a smooth 
rise and decay with comparable intensities in all eight 
detectors; this signature is inconsistent with a point 
source of photons. The second kind of event arises 
when the precipitation occurs at some distance from 
the spacecraft, so that it appears only in the detectors 
on one side of the satellite. Because the source of radi- 
ation is relatively nearby, however, the orbital motion 
of CGRO gives the event different profiles in differ- 
ent detectors. The third kind of event shows rapid 
variability with complex temporal structure, and can 
closely mimic a GRB. Such events can be recognized 
if they show characteristics from the first two classes, 
such as appearance in opposite-facing detectors or 
inconsistent time profiles between detectors. These 
events also exhibit a bremsstrahlung cutoff in their 
energy spectrum. 

If an event does not appear to be any of the above, 
we define background intervals by hand and estimate 
the source direction using our version of the BATSE 
LOCBURST software (see below). We also examine 
the profile of the event in the 4 DISCLA energy chan- 
nels. 

Events with directions estimated to be behind the 
Earth are classified as "Earth" events. We have not 
yet done a detailed analysis of this category, but it 
is likely to include electron precipitation events oc- 
curring below the spacecraft, events from other cate- 
gories for which we obtained poor direction estimates, 
and possibly terrestrial gamma ray flashes (TGRFs; 
Fishman et al. 1994a). 



Solar flares are identified based primarily on spec- 
tral softness and location. Typical solar flares have 
most of their counts in DISCLA channel 1, with fewer 
in channel 2, and fewer still in channel 3. In contrast, 
typical GRBs have most of their counts in channels 
2 and 3, with proportionally less signal in channel 1. 
Some spectrally hard solar flares are observed, how- 
ever, so events with directions consistent with the sun 



(within uncertainties) are classified as solar. For the 
search described here we did not compare laboratory 
trigger times with records of solar activity to separate 
hard solar flares from GRBs. 

Events which appear to be neither terrestrial nor 
solar and which have sufficient spectral hardness to be 
seen both in channels 2 and 3 are classified as GRB 
candidates. 

Events which do not make it into any of the previ- 
ous categories are classified as "unknown" , a category 
which includes all the low-energy (channel 1 only) 
events that are not obviously of terrestrial or solar 
origin. Events with significant counts in channels 1 
and 2 (25-100 keV) but not in channel 3 (100-300 
keV) are also included in this category. 

The classification of laboratory triggers is subjec- 
tive in cases where there is not an obvious indication 
of the nature of the event. To compare our classifica- 
tion methods with those of the BATSE team, we can 
use the 317 onboard triggers that were also laboratory 
triggers in our search of 345 days of data (see section 
H the remaining 256 triggers were not detected be- 
cause of gaps in the DISCLA data files). Table || 
shows how these onboard triggers were classified by 
the BATSE team and by us. 

Out of the 221 events that were identified as GRBs 
by the BATSE team, we classified 198 as GRBs, 10 as 
solar flares, 4 as magnetospheric events, and 9 in other 
categories (such as "unknown" or "Earth"). The 10 
events that were classified as GRBs by the BATSE 
team but as solar flares by us all had estimated direc- 
tions consistent with the sun. We classified them as 
solar flares because we could not argue that they were 
GRBs rather than hard solar flares. The 4 BATSE 
GRBs that we classified as magnetospheric events oc- 
curred while the spacecraft was in the vicinity of the 
SAA or at a maximum (or minimum) geographic lati- 
tude, where particle precipitation events are common. 
The remaining 9 BATSE GRBs that we put in other 
categories include very short (duration < 1 s) bursts 
that were difficult to classify due to low signal-to- 
noise, events for which we estimated the directions to 
be behind the Earth, and bursts for which there was 
an unusually low count rate in channel 3 (100-300 
keV). 

The results in Table |2| show that our classifica- 
tions agree with those of the BATSE team in most 
cases. When there is uncertainty between a GRB and 
a non-GRB origin, our classification is "conservative" 
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in that we tend towards the non-GRB classification. 
There were only 2 events out of 200 that we classified 
as GRBs but that the BATSE team classified as mag- 
netospheric particle precipitation events. Over the 
course of a year, our tendency to classify hard events 
occurring near the sun as solar flares introduces a bias 
against GRBs that occur in the plane of the ecliptic. 

For each of the solar flares, GRB candidates, and 
unknown events we estimate a source direction, in- 
tensity, and power-law spectral index using a modi- 
fied version of the BATSE LOCBURST code. This 
software uses the BATSE detector response matri- 
ces along with a model for the scattering of incident 
photons by the Earth's atmosphere to find the count 
rates expected from an event with intensity P (in pho- 
tons cm~ 2 s _1 above 10 keV), power-law spectral in- 
dex v, and source direction (9, <fi) in CGRO coordi - 
nates ( Fishman ct al. 1994bj ; Pendleton ct al. 1995 ). 
Let C^(P,v,6,(/)) denote the model count rates ex- 
pe cted in energy channel i of detector d, and let C\ 



der jote the measured (background subtracted) count 
rates with associated statistical measurement uncer- 
tainties &i t d- To estimate P, v, 9, and <fr for a given 
event, the software minimizes the following measure 
of goodness-of-fit, 



X 2 (P,i>,6,0)=Y, 

i . d 



ci-c* d (p,i>,i 



(3) 



The measured count rates used for this procedure are 
mean count rates obtained from the main (most in- 
tense) portion of the burst as selected by hand during 
our visual inspection of the laboratory triggers. 

The angular response of the BATSE detectors can 
allow for multiple, widely separated local minima in 
the x 2 parameter space, especially for weak bursts. 
For example, if a burst has most of its counts in just 
two detectors it can be equally consistent with two 
burst directions depending on the choice of the third 
most brightly illuminated detector. For weak bursts, 
the statistics may not be good enough to reliably dis- 
tinguish the third detector. Background variations 
(due to real source activity) in some detectors can 



also make it difficult to distinguish the third detec- 



tor. In such situations, \ 2 can be a strongly non- 
linear function of the observed rates. The errors on 
the estimated burst parameters cannot then be reli- 
ably estimated from the formal covariance matrix of 
the fit. 

To estimate errors on the model parameters for 



each burst, we produce 50 sets of "synthetic" burst 
rates obtained by drawing from a random distribution 
with the same means and variances as the measured 
rates. These 50 sets of synthetic count rates are then 
subjected to the same \ 2 minimization procedure as 
the real rates. The variances in the parameters ob- 
tained from the synthetic count rates are used to es- 
timate the uncertainty on the parameters (P, i>, 9, (j)) 
obtained from the measured count rates. 

Using the estimated mean intensity P we obtain a 
conversion factor from counts s _1 in the most brightly 
illuminated detector to units of ph cm~ 2 s _1 in the 
50-300 keV range. The peak flux and fluence of each 
event in physical units are then determined by mul- 
tiplying the corresponding measured counts by this 
conversion factor. The durations of events are char- 
acterized by the T50 and T90 duration measures, which 
are the time intervals during which the burst fluence 
increases from 25% to 75% and from 5% to 95% (re- 



spect i vely) of the total fluence ( Kouvcliotou ct al 
1993; Koshut et al. 1996). Uncertainties in the peak 



flux, fluence, and durations are derived from the un- 
certainties in P, z>, and the measured (background 
subtracted) count rates using the standard techniques 
for the propagation of small random errors (although 
the assumptions required by this method are not al- 
ways satisfied). 

We have attempted to ascertain how well the above 
procedures estimate burst intensities and directions 
by applying our methods to GRBs from the 3B cata- 
log ( Meegan et al. 1996| ). Because the burst intensi- 
ties and directions must be estimated simultaneously 
by folding a model through the detector and atmo- 
spheric response matrices, systematic errors in the 
inferred quantities can arise from the background sub- 
traction, the modeling of the event spectrum, spectral 
changes during the event, and detector calibration. 

To evaluate the accuracy of our intensity mea- 
surements, we applied our analysis procedure to 29 
GRBs from the 3B catalog ( [Vleegan ct al. 1996Q . We 
compared our peak fluxes (derived from the DISCLA 
data) to those obtained by the BATSE team (using 



the high resolution burst data types; see Pendleton et 
1996| ). In 13 bursts the two measurements agreed 



al 



to within the ler statistical uncertainties, and in 24 
bursts they agreed to within 2a. The largest dis- 
agreements (in terms of standard deviations) occur 
only in the most intense bursts, where the systematic 
errors are expected to dominate the statistical uncer- 
tainties. In those cases, the measurements disagree 
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by less than 30%. 

At least two sources of systematic errors are re- 
sponsible for the differences between our peak flux 
measurements (on intense bursts) and those obtained 
by Pendleton et al. (f996) for the 3B catalog. First, 
we model the incident photon spectrum as a power 
law, whereas Pendleton et al. (f996) allow for some 
curvature in the spectrum. Second, our peak flux 
measurement on the f024 ms time scale is based on 
the D1SCLA time bin that contains the most counts 
above background. The phase of this 1024 ms bin rel- 
ative to the "true" 1024 ms peak flux depends on the 
DISCLA sampling times; therefore, the "true" 1024 
ms peak flux could be spread over two DISCLA time 
bins. In contrast, Pendleton et al. (1996) use data 
with 64 ms time resolution to find which placement 
of a 1024 ms interval yields the highest peak flux on 
that time scale. 

A similar procedure was used to assess the accu- 
racy of our direction estimates. Comparison of a sam- 
ple of GRBs from the 3B catalog and hard solar flares 
which triggered onboard suggests that an additional 
systematic uncertainty of about 4° should be added 
in quadrature to our statistical direction uncertain- 
ties for events with emission between 50 and 300 keV. 
This is expected since our version of the LOCBURST 
code roughly corresponds (with minor improvements) 
to the version used to produce the BATSE IB cata- 
log (Fishman et al. 1994b). The improvements to 



LOCBURST that reduced systematic errors to 1.6° 
for the 3B catalog require more spectral information 
than is available in the DISCLA data and have not 
been incorporated into the analysis described in this 
paper. For the channel 1 only (25-50 keV) events, 
comparison with solar flares indicates that the addi- 
tional systematic uncertainty is about 6°. 

For faint events the background subtraction can be 
a substantial source of systematic error because it is 
not always clear what is background and what is low- 
level emission before or after the event. The duration 
estimates (and thus fluence estimates) are particu- 
larly sensitive to the choice of background intervals 
(see |Koshut et al. 1996| ). 



5. Results 

We have applied the search and analysis procedures 
described in sections || and [| to the DISCLA data 
taken between 1993 January 13 and 1993 December 
24. The corresponding range of Truncated Julian Day 



numbers is TJD 9001-9345 (TJD = Julian Day - 
2,440,000.5). Figure ^ shows a sky map which com- 
bines events from the GRB candidate (o), solar flare 
(*), and unknown (+) categories. The concentration 
of solar events in the ecliptic plane is clearly visible, 
as is a concentration of low-energy (+) events in the 
vicinity of Cyg X-l (£ = 71.3°, b = 3.0°). 

5.1. Gamma ray burst candidates 

Our search so far yields 91 non-triggered GRB can- 
didates. Tables || and || together are a catalog of these 
events. 

The first column of Table [| is a name which spec- 
ifies the approximate time of the event in the for- 
mat NTB yymmdd.ff, where yy is the year, mm is 
the month, and dd.ff is the day. The second column 
gives the time of the laboratory trigger expressed as 
the TJD and the seconds of day (SOD). The next 
three columns give the estimated source direction in 
equatorial (J2000) coordinates and its associated sta- 
tistical uncertainty. The full l.Ocr uncertainty in the 
direction estimates is obtained by combining the sta- 
tistical uncertainty and the 4° systematic uncertainty 
(see section Q) in quadrature. The next column gives 
the largest value (among the three time scales in Table 
|l|) of C max /C m ini which is the ratio of the maximum 
count rate achieved during the event to the minimum 
count rate required for detection in the 50-300 keV 
band. Events with C max /C m i n < 1.0 were detected in 
the "catch-all" search of the 25-300 keV band (search 
"c" of Table 0) and may represent spectrally softer 
GRBs that had too few counts to trigger in the 50- 
300 keV band. The next column gives the threshold 
number of counts in the 50-300 keV band, C m i n , and 
(in superscript) the search time scale which yielded 
the largest value of C max /C m i n . The next column 
lists the searches from Table [j] in which the event was 
detected. 

The last column of Table || gives the reasons in our 
estimation why the events did not trigger onboard the 
spacecraft. The notation "R" indicates the event oc- 
curred during the readout of a previous onboard trig- 
ger. "D" indicates that the event occurred while the 
onboard trigger was disabled due to passage through 
a region with a high probability of magnetospheric 
particle precipitation events. "F" indicates that the 
event was too faint to meet the onboard burst trig- 
ger threshold. "BB" indicates that the event failed to 
trigger onboard because it biased the onboard back- 
ground average. In one case, no experiment house- 
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keeping data were available so the state of the on- 
board trigger at the time of the event could not be 
determined. A "?" is entered for this event because 
we cannot determine whether the reason is "D" or 
"BB" . (See below for further discussion of these rea- 
sons and examples of events in each category.) 

Table |] gives the durations and intensities (in phys- 
ical units) of the GRB candidates. The first column 
gives the name of the GRB candidate. The next two 
columns give estimates of the T50 and T90 duration 
measures. Events with no entry in the T50 column 
had their Tgo duration estimated by eye. Uncertain- 
ties listed as 0.00 indicate that the uncertainty is less 
than the duration of one DISCLA time bin (1.024 s). 
The next column gives the peak flux in the 50-300 
keV range as measured on the 1024 ms time scale. 
The next column gives the 50-300 keV fluence esti- 
mate. 

Figure || shows a sky map of the direction estimates 
for the 91 untriggered GRB candidates. Events are 
shown as 1.0a error circles centered on the best- fit lo- 
cation. Using the sky exposure calculated in section 
3.1, the dipole and quadrupole moments of this direc- 
tion distribution in Galactic coordinates are (cos 6) = 
-0.001 ± 0.025 and (sin 2 b - 1/3) = -0.022 ± 0.021, 
where 8 is the angle between the burst direction 
and the Galactic center and b is the Galactic lon- 
gitude. These values are consistent with the values 
(cos 6) = 0.00 ± 0.06 and (sin 2 b - 1/3) = 0.00 ± 0.03 
expected from an isotropic distribution with the same 
number of bursts. 

The dipole and quadrupole moments (corrected 
for sky exposure) with respect to equatorial coor- 
dinates are (sini5) = 0.036 ± 0.027 and (sin S — 
1/3) = 0.074 ± 0.024. The dipole moment is con- 
sistent with that expected from an isotropic distri- 
bution, (sin S) — 0.00 ± 0.06. The quadrupole mo- 
ment appears to be only marginally consistent with 
the value (sin 2 S— 1/3) = 0.00 ± 0.03 expected for an 
isotropic distribution, indicating a weak concentra- 
tion of events in the direction of the celestial poles. 
This result may be due to our tendency to classify 
GRBs with directions consistent with the sun as hard 
solar flares (see section [|). 

The durations based on the Tgo interval of these 
events range from < 1.024 s to ~ 350 s. We have 
examined the non-triggered GRB candidates to see if 
any appear to be related to an onboard triggered GRB 
that occurred within 1 day of a non-triggered event. 
In a combined sample of 91 non-triggered GRB can- 



didates and 333 bursts from the 3B catalog (covering 
T JD 8995-9347) , we found 7 pairs of bursts occurring 
within 1 day of each other and having direction mea- 
surements compatible within la uncertainties. Only 
3 of these 7 pairs involved a non-triggered GRB can- 
didate. We do not consider this to be evidence that 
any of these pairs share the same burst source. We 
expect statistically to find 5-8 such pairs in a sam- 
ple of the same size drawn from bursts randomly dis- 
tributed uniformly in time and isotropically in space, 
with median location measurement errors of 8-10°. 
If one or more convincing pairs had been identified, 
however, they could have either been interpreted as 
burst repetition (see Wang 1994 ) or as bursts that 
extend knowledge of the Tgo distribution to longer 
durations. 

Figure ^| shows the integral peak flux distribution 
on the 1024 ms time scale for the untriggered GRB 
candidates. These events are concentrated at the faint 
end of the distribution, as expected for events which 
were generally too faint to cause an onboard trigger. 
When the non-triggered bursts are added to the trig- 
gered bursts for the same time period, the departure 
from the —3/2 power law slope expected from a ho- 
mogeneous distribution (in Euclidean space) remains 
evident. 

The top two panels of Figure [ll] show the differen- 
tial, n(P), and integral, N(> P), distributions of the 
peak fluxes, P, for the combined sample of 83 non- 
triggered GRB candidates and 233 onboard triggered 
GRBs that were detected on the 1024 ms time scale. 
The bottom panel shows the slope of the logarithmic 
number versus peak flux distribution, defined by 



s(P) 



d\ogN(> P) 
dlogP 



-P 



n(P) 
N{> P) ' 



(4) 



The dotted histograms in Figure |T(J show the distri- 
butions corrected for our laboratory trigger efficiency 
(see figure ||). The logarithmic slope s(P) is con- 
sistent with a value of —0.5 ± 0.1 at peak fluxes of 
0.15-0.35 ph cm" 2 s" 1 . 

Another measure of the inhomogeneity of the source 
distribution is reflected in the distribution of the 
V/Vmax statistic for these events, gi ven by (C max /C min 
flSchmidt, Higdon, fc Hcutcr 1988| ). For the 91 non- 
triggered GRB candidates, the average value (V/V max ) 
0.28 ± 0.03. This value is biased by the elimination 
of the strong (triggered) bursts; but it is of interest 
when considering how biases inherent in the onboard 
trigger mechanism could affect conclusions about the 
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spatial inhomogeneity of GRB sources. For exam- 
ple, if the onboard trigger's bias against slow-risers 
had significantly biased the value of (V/V max ) — 0.33 



± 0.01 obtained for the 3B catalog ( Meegan et al 
19S[6j) then the value obtained from the non-triggered 
GRB candidates alone could be expected to be much 
higher. For comparison, the value obtained using the 
cataloged bursts detected during the same time pe- 
riod is (V/V m ax} — 0.13 ± 0.02. For the combined 
non-triggered GRB candidates and cataloged bursts, 
(V/V max ) =0.18 ±0.020 

Figure ^ shows intensity profiles of some repre- 
sentative GRB candidates. Two adjacent plots are 
shown for each event. The plots on the left show 
the burst profile from the detector most brightly il- 
luminated by the burst. The plots on the right show 
the profile from the second most illuminated (or "sec- 
ond brightest") detector; they illustrate why some 
of the events did not cause an onboard burst trig- 
ger. The dashed lines show our estimate, of the 
background counts in each bin (fc) based on poly- 
nomial fits to data before and after the event. The 
dotted lines show the 5.5crs threshold level given by 

5.5^/ B^ . This level represents an "ideal" threshold 
estimate and it is in general different from the actual 
onboard background estimate that was in effect at the 
time the event occurred (see section ||) . The onboard 
background estimate has a statistical uncertainty re- 
sulting from the uncertainty in the mean count rate 
during the 17.408 s background accumulation inter- 
val. The dot-dashed line represents the threshold cor- 
responding to 5.5y B^ plus the uncertainty in the 
onboard trigger level arising from the onboard back- 
ground uncertainty. An event with peak counts just 
slightly above our "ideal" 5.5o\b threshold may fall be- 
low the onboard burst trigger threshold even though 
both thresholds are based on background estimates 
that are statistically consistent with each other. We 

classify events with peak counts less than 

plus the onboard threshold uncertainty as too faint 

to trigger onboard ( "F" in table ||) . 

The event in row (a) of Figure [□] failed to trigger 
onboard because it occurred during the read out pe- 
riod of a previous onboard burst trigger. The event 
in row (b) occurred while the onboard burst trigger 



was disabled due to spacecraft passage through a re- 
gion identified with a high probability of a false trigger 
due to atmospheric electron precipitation events. The 
events in rows (c) and (d) failed to trigger onboard 
because the onboard background estimate was biased 
upwards by slowly rising burst flux; these are exam- 
ples of slow- risers. Panels (e) and (f) show events that 
were too faint to meet the onboard trigger threshold. 

We estimate that the 91 GRB candidates failed 
to trigger onboard the spacecraft for the following 
reasons: 15 events occurred during the read out of 
a brighter event, 2 occurred while the onboard trig- 
ger was disabled for other reasons, 63 were below the 
5.5(7b threshold in the second brightest detector, and 
10 had a slow rise that modified the onboard back- 
ground estimate. One occurred during a time for 
which no spacecraft housekeeping data are available 
to determine the status of the onboard trigger. 

The onboard trigger mechanism's bias against slow- 
rising GRBs has been discussed by Lingenfelter and 
Higdon (1996). The 10 (possibly 11) events we find 
that failed to trigger onboard the spacecraft solely 
because of the slow-rising effect constitute 3.0% (pos- 
sibly 3.3%) of the total 332 GRBs that have been de- 
tected above the onboard threshold while the trigger 
was active. This is a lower fraction than estimated 
elsewhere ( Lingenfelter &: Higdon 1996| ). We note, 



3 These values are corrected from the erroneous values that ap- 
peared in earlier versions of this paper and in the version that 
appeared in ApJ, 491, 704 (1997). See erratum at end of paper. 



however, that our search algorithm is biased against 
faint events which rise on time scales longer than ~ 30 
s on the 1.024 s time scale (see Figure^). 

5.2. Unknown events 

The "unknown" category of laboratory triggers in- 
cludes all events which were not obviously of terres- 
trial or solar origin and which do not resemble a GRB. 
Most of the channel 1 only (low-energy, 25-50 keV) 
events fall into this category. The major problem with 
this class of events is that it is dominated by intensity 
fluctuations from Cyg X-l: of 799 events in the un- 
known category, 689 are consistent with this source 
(although they may not all be from Cyg X-l); see 
the clustering of events marked (+) in Figure [?]. If 
we remove all the events consistent with Cyg X-l we 
are left with the sky map shown in Figure ^|, where 
events are plotted as their l.Otr error circles. Although 
Figure |lj shows some general clustering toward the 
galactic center, there is no obvious clustering that 
would indicate the activity of any particular source. 

We find two events which can convincingly be at- 
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tributed to SGR 1806—20 based on intensity, spectral 
softness, and location. Both occur within one day of 
the onboard triggered emission from SGR 1806—20 
reported by Kouveliotou et al. (1994). Recent ac- 
tivity from this source suggests that more events 
from SGR 1806-20 (or other SGRs) may be de- 
tectable when this search is extended (Kouveliotou 
ct 41. 1996j ). 



6. Conclusions 

Our search of 345 days of archival BATSE data 
has uncovered a significant number of astronomically 
interesting transients. 

The 91 non-triggered GRB candidates detected (so 
far) by this search include some of the faintest GRBs 
ever observed. When combined with the bursts de- 
tected by the onboard trigger during the same 345 
days, these events extend knowledge of the peak flux 
distribution to values a factor of ~ 2 lower than 
the onboard detection threshold. Near the onboard 
trigger threshold, the combined sample is expected 
to be nearly complete (on the 1024 ms time scale). 
We find the logarithmic slope of the integral num- 
ber versus peak flux distribution to be —0.5 ±0.1 at 
peak fluxes of 0.15-0.35 ph cm ~ 2 s _1 after correct- 
ing for our laboratory trigger efficiency. The value of 
(V/V max ) = 0.18 ± 0.02 for the combined sample.^ 
We find no evidence for anisotropy in the direction 
distribution of these events. 

These results are consistent with those obtained 
from analyses of the published BATSE catalogs. The 
biases inherent in the onboard trigger mechanism do 
not appear to have significantly undermined its sam- 
pling of GRBs, at least for bursts with the character- 
istics our search can detect. 

The non-triggered GRB candidates add to the 
database of GRBs available for duration studies, searches 
for burst r epetition, and searches fo r gravitational 
lensing (see Fishman fc Meegan 1995 for an overview 
and references). The slow-risers and the more in- 
tense bursts (which occurred while the onboard trig- 
ger was disabled) will probably be the most useful 
non-triggered GRBs for such purposes. 

The low-energy (25-50 keV) events detected by our 
search arise from a variety of sources. While inten- 
sity fluctuations from Cyg X-l dominate this class 



of events, we find a significant number that must be 
due to other sources. Because of the difficulty in accu- 
rately estimating the source directions of these events, 
identification of the individual sources responsible for 
them depends on unique repetition patterns or tem- 
poral coincidences with other observations (as in the 
case of events from SGR 1806—20). The possibility 
remains that we may identify among these events new 
source activity or completely new burst sources. 

The effort to extend this search to cover the more 
than 5 years of remaining archival data is in progress. 

J. M. K. acknowledges support from a National Sci- 
ence Foundation Graduate Research Fellowship dur- 
ing the preliminary phase of this research and sub- 
sequent support from NASA Graduate Student Re- 
searchers Program Fellowship NGT8-52816. W. H. 
G. L. acknowledges support from NASA under grant 
NAG5-3804. C. K. acknowledges support from NASA 
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thank Scott Barthelmy and James Kuyper for sharing 
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4 This value is corrected from the erroneous value that appeared 
previously. See erratum at end of paper. 
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ERRATUM 



In the paper "A Search for Nontriggered Gamma- 
ray Bursts in the BATSE Data Base" by Kommers 
ct al. (ApJ, 491, 704, [1997]) the values of (V/V max ) 
that include onboard-triggered bursts are incorrect. 
Owing to a programming error, whenever a value of 
C m in was available for a burst listed in the 4B catalog, 
that catalog value (appropriate for the onboard burst 
trigger) overwrote the value appropriate for our more 
sensitive off-line search. Thus the values of V/V max 
were overestimated by a factor of <~ 2 3 / 2 for the bursts 
that were triggered onboard. This error does not af- 
fect the values of C max / C m i n and C m i n for the non- 
triggered GRBs listed in Table 3, nor does it affect 
the value of {V/V max } = 0.28±0.03 computed for the 
non-triggered (only) GRB sample. 

The correct value of (V/V max ) for the onboard- 
triggered bursts detected by our off-line search is 
(V/V max ) = 0.13 ± 0.02. For the combined non- 
triggered GRB candidates and onboard-triggered bursts, 
{V/Vmax) = 0.18 ± 0.02. This value is significantly 
lower than the value of {V/V max } — 0.33 ± 0.01 ob- 
tained for the 3B catalog. These corrections, there- 
fore, significantly strengthen the conclusion of the pa- 
per, that the BATSE onboard trigger has not missed 
a significant population of faint gamma-ray bursts ow- 
ing to the biases inherent in the onboard burst trigger 
mechanism. 
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Fig. 1. — Schematic diagram of the laboratory (off- 
line search) trigger criteria. If the significances of the 
fluctuations measured in the two most brightly illu- 
minated detectors lie in the non-shaded area A, our 
search will flag that time bin as a laboratory trigger. 

Fig. 2. — Integral distribution of Sd(k) for the data of 
May 27, 1993 (TJD 9135) in detector 0. The dotted 
line shows the expectation for a Gaussian distribution 
with zero mean and unit variance. The excess counts 
are attributed both to deficiencies in the background 
estimate and to the activity of real sources. 

Fig. 3. — Off-line trigger efficiency. The solid line 
shows the efficiency of our off-line search algorithm 
for detecting an event with a given peak flux in the 
50-300 keV range on the 1024 ms time scale. The 
long-dashed line shows the trigger efficiency for our 
off-line search on the 4096 ms time scale. The short- 
dashed line shows the trigger efficiency from the IB 
catalog. 

Fig. 4. — Sky exposure for our 50-300 keV off-line 
trigger search as a function of declination, assuming 
a search spanning 345 days. 

Fig. 5. — Simplified profile of a "slow-riser". It 
reaches maximum significance s by rising with a slope 
r given in sigmas per second. 

Fig. 6. — Contours showing the probability of de- 
tecting a "slow-riser" with a given peak significance 
(relative to the "true" background) and a given slope. 
Thick contours apply to our off-line search algorithm, 
and thin contours apply to the onboard burst trigger. 

Fig. 7. — Sky map of all non-terrestrial non-triggered 
events in Galactic coordinates. The GRB candidate 
(o), solar flare (*), and unknown (+) categories are 
shown. The concentration of solar events in the eclip- 
tic plane is clearly visible, as is a concentration of 
low-energy (+) events in the vicinity of Cyg X-l 
(£ = 71°.3, b = 3°.0). 

Fig. 8. — Sky map of 91 non-triggered GRB candi- 
dates, shown as la error circles in Galactic coordi- 
nates. 



Fig. 9. — Integral number versus peak flux distri- 
bution of 91 non-triggered GRB candidates (dotted 
line). No corrections for trigger efficiency have been 
applied. The distribution for GRBs from the 3B cat- 
alog detected during the same time period is also 
shown (dashed line), as is that for the combined sam- 
ple (solid line). 



Fig. 10. — Peak flux distributions for 83 non-triggered 
GRB candidates combined with the 233 onboard trig- 
gered events (1024 ms time scale) from the same time 
period. The solid histogram shows the observed num- 
bers and the dotted histogram shows the numbers 
corrected for laboratory trigger efficiency. n(P) is the 
differential distribution, N(P) is the integral distribu- 
tion, and s(P) is the slope of the logarithmic number 
versus peak flux distribution. 



Fig. 11. — Intensity profiles of selected GRB candi- 
dates. Two panels are used for each event. Those on 
the left represent the count rates observed in the de- 
tector most brightly illuminated by the burst. Those 
on the right represent the count rates in the second 
most illuminated detector. Dashed lines represent our 
background estimate derived from polynomial fits to 
data before and after the event. Dotted lines repre- 
sent an "ideal" onboard trigger threshold based on 
our background estimates; the actual onboard trigger 
threshold is in general different and is based on an 
average count rate that is recomputed every 17.408 s. 
The dot-dashed line represents the "ideal" onboard 
threshold plus the uncertainty arising from the statis- 
tical uncertainty in onboard background average. The 
burst in row (a) occurred during the read out period 
of a more intense event, and that in row (b) occurred 
while the onboard trigger was disabled. The bursts 
in rows (c) and (d) modified the onboard background 
average. Rows (e) and (f) show events that were too 
faint cause an onboard burst trigger. 



Fig. 12. — Low-energy (25-50 keV) events plotted as 
la error circles in Galactic coordinates (those consis- 
tent with Cyg X-l excluded). 
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Table 1 

Parameters of the time series formed from the DISCLA data (see text). 



Search 


Time Bin 
Duration (s) 


Energy 
Channels 


M 


£ 


N b 


Statistical Detections 
per day 


a 


1.024 


1 


2.5 


4.0 


20 


0.021 


b 


1.024 


2+3 


2.5 


4.0 


20 


0.021 


c 


1.024 


1+2+3 


2.5 


4.0 


20 


0.021 


d 


4.096 


1 


2.5 


4.0 


15 


0.005 


e 


4.096 


2+3 


2.5 


4.0 


15 


0.005 


f 


8.192 


1 


2.5 


4.0 


5 


0.002 


g 


8.192 


2+3 


2.5 


4.0 


5 


0.002 



16 



Table 2 

Classification matrix of 317 onboard triggers that were also detected by our off-line search, 
rows indicate the classification assigned by the batse team, and columns indicate the 

classification assigned by us (see text) . 





Off-line 


Off-line 


Off-line 


Off-line 




GRB Candidate 


Solar Flare 


Magnetospheric 


Other 


BATSE GRB 


198 


10 


4 


9 


BATSE Solar Flare 





50 


3 





BATSE Magnetospheric 


2 


2 


30 


4 


BATSE Unclassified 











5 
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Table 3 

Times, source directions, C max /C min , searches triggered, and reasons for eluding the onboard 

TRIGGER FOR 91 NON-TRIGGERED GRB CANDIDATES. 



Name 


Time 

1 J L).S) 


RA 

\ ) 


Dec 
I ) 


Err. 

( ) 


Cmax 

G min 


C min 


Searches 
Triggered 


Reason 
Non- triggered 


NTB 930118.74 


9005:64425.6 


219.3 


-32.9 


0.5 


31.0 


349 4 


abedefg 


R 


NTB 930211.88 


9029:76428.9 


285.5 


20.0 


7.7 


3.5 


425 8 


bedefg 


F 


NTB 930216.63 


9034:54956.7 


42.5 


-8.7 


3.6 


2.6 


447 8 


bedeg 


F 


NTB 930217.80 


9035:69525.6 


311.7 


-12.0 


28.2 


0.9 


286 4 


c 


F 


NTB 930225.86 


9043:75141.2 


336.6 


35.2 


38.4 


1.1 


155 1 


abc 


F 


NTB 930227.83 


9045:71910.5 


241.2 


81.0 


11.0 


2.1 


377 8 


cog 


F 


NTB 930228.85 


9046:73728.6 


18.6 


2.5 


8.3 


2.2 


477 8 


g 


F 


NTB 930302.20 


9048:17613.4 


214.1 


42.5 


6.2 


2.4 


483 8 


g 


F 


NTB 930303.65 


9049:56726.7 


0.0 


-50.8 


3.3 


3.2 


501 8 


bceg 


F 


NTB 930305.70 


9051:60923.0 


275.2 


59.8 


3.5 


3.2 


330 8 


bceg 


F 


NTB 930307.54 


9053:46677.1 


94.2 


-15.5 


12.7 


1.7 


411 8 


be 


F 


NTB 930308.30 


9054:26710.7 


88.7 


-31.6 


3.5 


3.8 


399 8 


bceg 


F 


NTB 930310.08 


9056: 7334.5 


333.5 


-57.2 


5.0 


2.6 


321 8 


abedefg 


F 


NTB 930315.46 


9061:40070.8 


232.5 


-32.1 


3.9 


2.0 


395 8 


g 


F 


NTB 930316.74 


9062:64295.0 


314.1 


-87.1 


19.8 


1.6 


132 1 


be 


F 


NTB 930318.18 


9064:15764.1 


16.2 


44.0 


8.9 


2.4 


425 8 


bceg 


F 


NTB 930320.94 


9066:81767.0 


228.3 


72.7 


15.5 


7.2 


436 s 


bcefg 


R 


NTB 930325.65 


9071:56254.1 


46.9 


46.8 


5.8 


2.8 


383 8 


bceg 


F 


NTB 930327.46 


9073:40594.6 


178.7 


-5.5 


13.2 


3.4 


553 8 


be 


F 


NTB 930330.91 


9076:79132.8 


109.2 


-54.4 


2.0 


1.7 


481 8 


cfg 


F 


NTB 930403.84 


9080:73239.2 


245.8 


-59.6 


6.2 


3.0 


351 8 


abedefg 


BB 


NTB 930409.13 


9086:11442.8 


315.6 


68.8 


4.0 


3.3 


334 8 


abedef 


BB 


NTB 930409.91 


9086:78639.7 


275.4 


-17.7 


0.9 


12.8 


501 8 


abedefg 


R 


NTB 930410.76 


9087:65711.8 


182.3 


-45.3 


7.3 


3.0 


380 8 


abeefg 


F 


NTB 930416.56 


9093:48602.2 


115.6 


-21.0 


6.2 


3.1 


389 8 


abedefg 


F 


NTB 930417.78 


9094:68020.4 


190.5 


8.3 


6.7 


1.7 


351 8 


abedf 


F 


NTB 930421.11 


9098:10164.9 


23.7 


19.3 


7.0 


2.0 


520 8 


g 


F 


NTB 930422.58 


9099:50820.2 


50.4 


-7.5 


2.1 


4.9 


315 8 


be 


? 


NTB 930424.45 


9101:38903.4 


230.3 


-55.7 


0.4 


21.2 


523 8 


abede 


D 


NTB 930424.97 


9101:84156.0 


254.1 


68.6 


3.2 


3.5 


332 8 


abeefg 


BB 


NTB 930426.48 


9103:41832.1 


33.2 


-81.9 


5.4 


4.4 


133 1 


abedefg 


R 


NTB 930427.59 


9104:51155.6 


60.1 


34.3 


11.5 


2.2 


453 8 


abd 


F 


NTB 930429.75 


9106:65094.3 


35.7 


-25.8 


5.4 


5.7 


344 8 


bedefg 


BB 


NTB 930501.34 


9108:29834.4 


181.4 


-32.9 


31.7 


1.1 


114 1 


ab 


F 


NTB 930506.63 


9113:55244.4 


259.8 


35.4 


8.4 


2.4 


548 8 


eg 


F 


NTB 930508.95 


9115:82814.6 


82.9 


41.4 


7.6 


3.4 


494 8 


bee 


D 


NTB 930513.98 


9120:85533.8 


169.0 


12.0 


39.5 


1.5 


384 8 


be 


F 


NTB 930519.39 


9126:34288.8 


272.4 


-24.2 


83.8 


1.4 


134 1 


be 


F 


NTB 930612.63 


9150:55165.1 


254.5 


34.3 


2.9 


12.1 


277 4 


abede 


R 


NTB 930616.27 


9154:23806.6 


179.5 


1.8 


4.4 


11.7 


165 1 


abede 


R 


NTB 930617.23 


9155:20027.0 


244.9 


-12.7 


6.1 


2.5 


299 8 


bceg 


F 


NTB 930626.94 


9164:81935.0 


342.6 


-34.9 


3.3 


3.6 


512 8 


bedefg 


BB 


NTB 930630.71 


9168:61420.2 


70.0 


38.1 


6.9 


3.0 


438 8 


bedeg 


F 


NTB 930701.62 


9169:54302.9 


226.4 


39.4 


2.5 


4.2 


338 8 


abede 


F 


NTB 930705.64 


9173:55983.3 


195.3 


-59.5 


7.2 


2.0 


334 8 


bcefg 


F 


NTB 930717.20 


9185:18101.4 


184.3 


57.8 


10.3 


1.5 


364 4 


bed 


F 


NTB 930717.98 


9185:85357.7 


200.5 


-66.4 


3.5 


3.3 


374 8 


abc 


F 


NTB 930722.84 


9190:73297.6 


310.5 


-48.2 


3.5 


2.1 


412 8 


eg 


F 


NTB 930728.54 


9196:47072.9 


90.8 


19.0 


21.0 


3.3 


510 8 


d 


R 


NTB 930804.71 


9203:61858.5 


32.7 


66.6 


5.0 


3.7 


317 8 


bcefg 


F 


NTB 930811.62 


9210:53728.4 


347.5 


65.0 


3.6 


2.9 


325 8 


be 


F 


NTB 930812.27 


9211:23904.9 


198.0 


-27.8 


2.4 


4.9 


409 8 


bedef 


BB 


NTB 930813.76 


9212:65850.5 


76.9 


77.1 


8.1 


3.7 


336 8 


abedefg 


F 


NTB 930816.67 


9215:58569.4 


155.1 


53.8 


5.9 


3.5 


429 s 


bceg 


F 


NTB 930820.76 


9219:65885.3 


62.7 


36.4 


7.6 


2.8 


386 s 


bceg 


F 


NTB 930821.64 


9220:56096.9 


148.0 


-40.4 


11.7 


1.2 


123 1 


be 


F 


NTB 930825.48 


9224:41775.3 


59.7 


63.1 


21.4 


1.5 


298 4 


eg 


F 


NTB 930827.60 


9226:51963.0 


349.4 


68.5 


19.0 


2.5 


506 8 


bceg 


F 


NTB 930902.45 


9232:39001.8 


224.2 


20.8 


6.6 


2.1 


364 8 


acg 


F 


NTB 930918.46 


9248:39913.1 


275.3 


-81.0 


11.8 


2.6 


310 8 


bceg 


F 


NTB 930921.84 


9251:73133.2 


75.0 


-35.7 


25.2 


1.3 


292 8 


bee 


F 


NTB 930924.37 


9254:32251.0 


98.1 


-9.6 


21.5 


3.3 


403 8 


bede 


F 


NTB 930928.93 


9258:81199.3 


260.5 


-67.6 




2.2 


529 8 


bef 


F 


NTB 930928.94 


9258:81393.3 


79.8 


45.2 


c 


1.7 


329 8 


f 


F 


NTB 931001.06 


9261: 5859.5 


207.2 


12.6 


1.4 


5.6 


582 8 


abede 


R 


NTB 931001.72 


9261:62917.3 


8.5 


12.5 


2.2 


3.4 


473 8 


bcefg 


BB 


NTB 931007.20 


9267:17771.7 


311.3 


-5.2 


23.7 


2.2 


368 8 


bg 


F 


NTB 931007.33 


9267:29319.3 


355.1 


36.4 


48.7 


1.7 


363 8 


bg 


F 



Table 3 — Continued 



Name 


Time 
(TJD:s) 


RA 

(°) 


Dec 

(°) 


Err. 

O 


<?min 


C min 


Searches 
Triggered 


Reason 
Non-triggered 


NTB 931008.05 


9268: 5024.9 


26.4 


-65.6 


30.5 


1.2 


190 1 


be 


F 


NTB 931008.63 


9268:54916.8 


42.3 


42.8 


3.7 


4.7 


492 8 


ccg 


BB 


NTB 931011.96 


9271:83718.3 


248.4 


63.9 


23.4 


1.2 


373 4 


be 


F 


NTB 931014.08 


9274: 7552.2 


65.5 


72.0 


55.6 


1.5 


148 1 


be 


F 


NTB 931017.22 


9277:19221.7 


185.7 


-67.5 


18.3 


1.8 


448 8 


b 


F 


NTB 931020.10 


9280: 8697.5 


285.0 


16.0 


14.0 


1.7 


447 8 


c 


F 


NTB 931025.93 


9285:80462.5 


173.2 


15.9 


4.8 


2.7 


438 8 


bceg 


F 


NTB 931031.23 


9291:20519.6 


173.2 


63.4 


17.2 


3.0 


163 1 


bee 


R 


NTB 931106.48 


9297:42228.4 


34.3 


69.7 


3.7 


3.7 


406 8 


bcefg 


BB 


NTB 931106.90 


9297:78310.6 


185.7 


-38.0 


16.1 


2.7 


279 4 


be 


R 


NTB 931107.31 


9298:26896.0 


106.0 


-27.7 


12.6 


2.9 


578 8 


abede 


BB 


NTB 931111.71 


9302:61941.4 


244.6 


49.5 


0.8 


9.2 


352 s 


abedefg 


R 


NTB 931113.04 


9304: 3669.7 


50.9 


-40.2 


10.1 


2.1 


383 8 


bee 


F 


NTB 931115.77 


9306:66557.6 


279.1 


27.9 


11.0 


3.1 


401 4 


abecg 


R 


NTB 931125.86 


9316:74847.4 


177.7 


-77.4 


2.4 


12.0 


529 8 


abedefg 


R 


NTB 931206.45 


9327:39631.5 


174.1 


-14.7 


15.0 


3.1 


404 8 


bceg 


F 


NTB 931209.89 


9330:77266.1 


20.0 


-37.4 


3.0 


3.3 


443 8 


bceg 


F 


NTB 931215.12 


9336:10491.5 


326.0 


-3.5 


12.2 


2.3 


298 4 


abd 


F 


NTB 931220.16 


9341:13826.7 


120.2 


45.7 


8.5 


3.0 


414 8 


bee 


F 


NTB 931220.73 


9341:63345.3 


50.0 


22.0 


3.9 


7.3 


435 8 


abedefg 


R 


NTB 931222.11 


9343:10361.5 


278.9 


29.5 


5.4 


1.9 


565 s 


g 


F 


NTB 931222.82 


9343:70972.1 


209.6 


-35.5 


14.8 


6.2 


178 1 


abecg 


R 


NTB 931223.07 


9344: 6589.6 


212.8 


40.0 


23.2 


1.5 


200 1 


be 


F 
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Table 4 

Durations and intensities of non-triggered GRB candidates. 



Name 


^ 50 


J 90 




r lucncc 






(ph cm 2 s 1 ) 


(ergs cm ) 


NTB 930118.74 


2.05 ± 0.00 


4.10 ± 0.00 


4.53 ± 0.14 


(3.00 ± 0.19)E-6 


NTB 930211.88 


19.46 ± 2.90 


58.37 ± 6.23 


0.26 ± 0.05 


(8.46 ± 1.95)E-7 


NTB 930216.63 


23.55 ± 0.00 


40.96 ± 8.44 


0.28 ± 0.07 


(6.34 ± 1.54)E-7 


NTB 930217.80 




3.07 ± 1.02 


0.18 ± 0.10 


(5.86 ± 3.50)E-8 


NTB 930225.86 




4.10 ± 1.02 


0.18 ± 0.07 


(7.72 ± 2.94)E-8 


NTB 930227.83 


37.89 ± 2.29 


57.34 ± 7.24 


0.17 ± 0.05 


(5.46 ± 1.36)E-7 


NTB 930228.85 


26.62 ± 3.24 


91.14 ± 16.67 


0.16 ± 0.04 


(9.25 ± 1.98)E-7 


NTB 930302.20 


21.50 ± 3.24 


57.34 ± 15.50 


0.26 ± 0.09 


(1.44 ± 0.44)E-6 


NTB 930303.65 


32.77 ± 1.45 


70.66 ± 7.17 


0.20 ± 0.05 


(7.64 ± 1.42)E-7 


NTB 930305.70 


5.12 ± 1.45 


12.29 ± 1.45 


0.59 ± 0.19 


(8.88 ± 3.73)E-7 


NTB 930307.54 


5.12 ± 1.02 


22.53 ± 16.79 


0.18 ± 0.19 


(2.26 ± 2.35)E-7 


NTB 930308.30 


18.43 ± 1.45 


38.91 ± 2.29 


0.43 ± 0.11 


(2.18 ± 0.53)E-6 


NTB 930310.08 


10.24 ± 0.00 


34.82 ± 23.57 


0.35 ± 0.05 


(4.16 ± 0.68)E-7 


NTB 930315.46 


15.36 ± 2.29 


35.84 ± 4.22 


0.23 ± 0.05 


(7.61 ± 1.81)E-7 


NTB 930316.74 




3.07 ± 1.02 


0.22 ± 0.07 


(8.05 ± 2.47)E-8 


NTB 930318.18 


15.36 ± 2.05 


35.84 ± 6.87 


0.17 ± 0.04 


(3.97 ± 0.94)E-7 


NTB 930320.94 


25.60 ± 0.00 


59.39 ± 3.07 


0.96 ± 0.15 


(4.09 ± 0.62)E-6 


NTB 930325.65 


10.24 ± 1.45 


20.48 ± 5.12 


0.24 ± 0.06 


(5.13 ± 1.42)E-7 


NTB 930327.46 


12.29 ± 3.69 


46.08 ± 31.89 


0.29 ± 0.07 


(7.86 ± 2.88)E-7 


NTB 930330.91 


12.29 ± 2.29 


52.22 ± 14.37 


0.24 ± 0.06 


(4.99 ± 1.47)E-7 


NTB 930403.84 


13.31 ± 0.00 


92.16 ± 10.69 


0.26 ± 0.04 


(6.96 ± 0.93)E-7 


NTB 930409.13 


14.34 ± 1.02 


39.94 ± 9.27 


0.40 ± 0.05 


(7.91 ± 1.20)E-7 


NTB 930409.91 


12.29 ± 1.02 


53.25 ± 0.00 


1.19 ± 0.06 


(3.06 ± 0.34)E-6 


NTB 930410.76 


24.58 ± 1.02 


74.75 ± 3.07 


0.29 ± 0.05 


(1.09 ± 0.16)E-6 


NTB 930416.56 




20.48 ± 5.12 


0.21 ± 0.04 


(3.54 ± 0.58)E-7 


NTB 930417.78 


5.12 ± 1.45 


35.84 ± 26.62 


0.21 ± 0.06 


(2.05 ± 0.82)E-7 


NTB 930421.11 


47.10 ± 3.24 


(1.29 ± 0.23)E+2 


0.15 ± 0.04 


(1.28 ± 0.26)E-6 


NTB 930422.58 


15.36 ± 1.02 


44.03 ± 17.41 


0.27 ± 0.03 


(1.20 ± 0.13)E-6 


NTB 930424.45 


41.98 ± 1.02 


(1.06 ± 0.03)E+2 


1.28 ± 0.06 


(1.38 ± 0.08)E-5 


NTB 930424.97 


12.29 ± 1.02 


32.77 ± 15.36 


0.33 ± 0.04 


(8.77 ± 1.30)E-7 


NTB 930426.48 




3.07 ± 1.02 


0.70 ± 0.10 


(1.66 ± 0.27)E-7 


NTB 930427.59 


21.50 ± 1.45 


46.08 ± 7.17 


0.19 ± 0.06 


(4.45 ± 1.30)E-7 


NTB 930429.75 


43.01 ± 1.45 


(1.11 ± 0.15)E+2 


0.38 ± 0.04 


(2.24 ± 0.25)E-6 


NTB 930501.34 




2.05 ± 1.02 


0.23 ± 0.14 


(5.71 ± 3.59)E-8 


NTB 930506.63 


34.82 ± 8.26 


82.94 ± 5.22 


0.18 ± 0.05 


(9.14 ± 2.73)E-7 


NTB 930508.95 


5.12 ± 1.45 


22.53 ± 15.39 


0.24 ± 0.05 


(4.17 ± 1.39)E-7 


NTB 930513.98 




41.98 ± 10.24 


0.17 ± 0.05 


(3.67 ± 1.07)E-7 


NTB 930519.39 




1.02 ± 0.00 


0.35 ± 0.29 


(1.15 ± 0.96)E-7 


NTB 930612.63 


6.14 ± 2.05 


50.18 ± 9.22 


1.33 ± 0.08 


(1.23 ± 0.42)E-6 


NTB 930616.27 




1.02 ± 0.00 


3.21 ± 0.77 


(7.97 ± 1.94)E-7 


NTB 930617.23 


51.20 ± 1.02 


73.73 ± 5.12 


0.29 ± 0.06 


(1.35 ± 0.28)E-6 


NTB 930626.94 


14.34 ± 1.02 


61.44 ± 38.01 


0.27 ± 0.05 


(8.29 ± 1.25)E-7 


NTB 930630.71 




18.43 ± 5.12 


0.23 ± 0.05 


(3.86 ± 0.71)E-7 


NTB 930701.62 


20.48 ± 1.02 


35.84 ± 3.07 


0.45 ± 0.05 


(1.07 ± 0.14)E-6 


NTB 930705.64 


12.29 ± 1.02 


48.13 ± 16.38 


0.20 ± 0.05 


(4.10 ± 0.95)E-7 


NTB 930717.20 


2.05 ± 1.45 


3.07 ± 1.02 


0.21 ± 0.09 


(1.37 ± l.ll)E-7 


NTB 930717.98 


6.14 ± 1.45 


17.41 ± 3.07 


0.32 ± 0.05 


(3.93 ± 1.09)E-7 


NTB 930722.84 


26.62 ± 2.29 


72.70 ± 9.44 


0.20 ± 0.05 


(8.37 ± 1.62)E-7 


NTB 930728.54 


66.56 ± 15.22 


(3.43 ± 0.21)E+2 


0.29 ± 0.05 


(3.60 ± 0.92)E-6 


NTB 930804.71 


19.46 ± 1.02 


49.15 ± 16.42 


0.25 ± 0.04 


(9.05 ± 1.10)E-7 


NTB 930811.62 


18.43 ± 2.29 


45.06 ± 4.22 


0.29 ± 0.06 


(7.27 ± 1.71)E-7 


NTB 930812.27 


28.67 ± 2.29 


77.82 ± 3.24 


0.39 ± 0.04 


(1.97 ± 0.25)E-6 


NTB 930813.76 


(1.02 ± 0.02)E+2 


(1.36 ± 0.02)E+2 


0.29 ± 0.05 


(1.16 ± 0.19)E-6 


NTB 930816.67 


8.19 ± 1.02 


24.58 ± 4.10 


0.33 ± 0.10 


(6.95 ± 2.21)E-7 


NTB 930820.76 


7.17 ± 1.02 


45.06 ± 23.55 


0.22 ± 0.05 


(4.03 ± 0.98)E-7 


NTB 930821.64 




2.05 ± 1.02 


0.30 ± 0.10 


(9.30 ± 3.26)E-8 


NTB 930825.48 


4.10 ± 2.29 


20.48 ± 5.12 


0.14 ± 0.06 


(1.58 ± 1.07)E-7 


NTB 930827.60 


15.36 ± 3.07 


48.13 ± 15.36 


0.22 ± 0.07 


(5.23 ± 1.80)E-7 


NTB 930902.45 


28.67 ± 10.69 


74.75 ± 16.51 


0.16 ± 0.04 


(7.53 ± 3.05)E-7 


NTB 930918.46 


16.38 ± 4.10 


54.27 ± 14.34 


0.20 ± 0.05 


(5.54 ± 1.95)E-7 


NTB 930921.84 


12.29 ± 3.24 


56.32 ± 18.43 


0.32 ± 0.12 


(3.12 ± 1.47)E-7 


NTB 930924.37 


8.19 ± 1.02 


28.67 ± 17.17 


0.22 ± 0.05 


(5.58 ± 1.23)E-7 


NTB 930928.93 


2.05 ± 1.02 


4.10 ± 2.05 


0.19 ± 0.06 


(6.92 ± 4.24)E-8 


NTB 930928.94 


14.34 ± 1.45 


28.67 ± 2.29 


0.13 ± 0.04 


(3.75 ± 1.04)E-7 


NTB 931001.06 


4.10 ± 0.00 


12.29 ± 3.07 


0.74 ± 0.08 


(8.66 ± 0.93)E-7 


NTB 931001.72 


16.38 ± 2.29 


82.94 ± 8.44 


0.28 ± 0.06 


(9.18 ± 1.98)E-7 


NTB 931007.20 


5.12 ± 1.02 


35.84 ± 5.12 


0.23 ± 0.07 


(3.62 ± 1.27)E-7 


NTB 931007.33 




8.19 ± 4.10 


0.43 ± 0.74 


(4.94 ± 8.59)E-7 
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Table 4 — Continued 



T 50 T 90 Peak Flux Fluonco 

(s) (s) (ph cm~ 2 s~ 1 ) (ergs cm -2 ) 



NTB 931008.05 








2.05 ± 1.02 


0.28 


± 





11 


(1.19 


± 


0.46)E-7 


NTB 931008.63 


66.56 


± 


1.45 


(1.98 ± 0.28)E+2 


0.26 


± 





04 


(3.34 


± 


0.29)E-6 


NTB 931011.96 


11.26 


± 


3.24 


50.18 ± 18.55 


0.20 


± 





11 


(2.00 


± 


1.27)E-7 


NTB 931014.08 








4.10 ± 2.05 


0.36 


± 





14 


(1.24 


± 


0.52)E-7 


NTB 931017.22 


9.22 


± 


5.22 


41.98 ± 23.57 


0.16 


± 





06 


(2.72 


± 


1.82)E-7 


NTB 931020.10 


21.50 


± 


3.24 


73.73 ± 2.29 


0.14 


± 





05 


(4.52 


± 


1.52)E-7 


NTB 931025.93 


20.48 


± 


9.22 


49.15 ± 15.66 


0.21 


± 





05 


(5.15 


± 


2.55)E-7 


NTB 931031.23 








1.02 ± 0.00 


0.40 


± 





11 


(8.95 


± 


2.43)E-8 


NTB 931106.48 


31.74 


± 


2.29 


90.11 ± 1.45 


0.22 


± 





04 


(1.41 


± 


0.21)E-6 


NTB 931106.90 








11.26 ± 4.10 


0.32 


± 





06 


(2.66 


± 


0.50)E-7 


NTB 931107.31 


30.72 


± 


3.24 


55.30 ± 9.27 


0.41 


± 





11 


(1.57 


± 


0.40)E-6 


NTB 931111.71 


58.37 


± 


3.07 


(1.86 ± 0.01)E+2 


0.77 


± 





04 


(7.24 


± 


0.59)E-6 


NTB 931113.04 


12.29 


± 


1.02 


52.22 ± 16.67 


0.21 


± 





05 


(4.15 


± 


0.93)E-7 


NTB 931115.77 


14.34 


± 


1.45 


26.62 ± 2.29 


0.32 


± 





06 


(3.15 


± 


0.76)E-7 


NTB 931125.86 


5.12 


± 


1.02 


14.34 ± 1.45 


0.83 


± 





05 


(1.27 


± 


0.27)E-6 


NTB 931206.45 


11.26 


± 


1.02 


29.70 ± 8.44 


0.26 


± 





05 


(5.43 


± 


l.ll)E-7 


NTB 931209.89 


34.82 


± 


1.02 


90.11 ± 4.58 


0.27 


± 





06 


(1.44 


± 


0.21)E-6 


NTB 931215.12 


3.07 


± 


2.05 


13.31 ± 11.31 


0.28 


± 





07 


(2.25 


± 


1.61)E-7 


NTB 931220.16 


4.10 


± 


1.45 


11.26 ± 4.10 


0.49 


± 





11 


(6.84 


± 


2.81)E-7 


NTB 931220.73 


8.19 


± 


0.00 


33.79 ± 4.58 


0.56 


± 





05 


(9.78 


± 


0.89)E-7 


NTB 931222.11 


18.43 


± 


3.24 


58.37 ± 12.83 


0.18 


± 





06 


(6.86 


± 


1.93)E-7 


NTB 931222.82 








2.05 ± 1.02 


1.75 


± 





24 


(4.40 


± 


0.65)E-7 


NTB 931223.07 








3.07 ± 1.02 


0.26 


± 





19 


(9.55 


± 


7.28)E-8 



21 





Fig. 2.— Fig. 4.— 
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Fig. 9.— 
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